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A  COMPARISON  OF  QUANTUM,  CLASSICAL,  AND  SEMICLASSICAL 
DESCRIPTIONS  OF  A  MODEL,  COLL  INEAR,  INELASTIC  COLLISION 
OF  TWO  DIATOMIC  MOLECULES 

I.  Introduction 

Molecular  dynamics  simulations  involving  a  relatively  large  number  of 
molecules  are  becoming  an  increasingly  powerful  tool  for  the  study  of  both 
gas  and  condensed  phase  phenomena.  There  has  been  considerable  progress  in 
the  ability  to  simulate  a  variety  of  equilibrium  and  nonequilibrium  situa¬ 
tions  and  in  calculating,  for  example,  transport  properties  (I).  This 
progress  is  due  both  to  the  improved  computational  facilities  and  improved 
methods  of  doing  the  calculations. 

There  is  a  now  interest  in  using  molecular  dynamics  to  describe  the 
interaction  of  a  shock  wave  and  a  condensed  phase  material.  The  motivation 
for  this  approach  is  based  on  the  observation  that  the  region  immediately 
behind  the  shock  wave  is  not  in  equilibrium.  Thus  molecular  dynamics  might 
be  used  when  the  usual  fluid  descriptions,  which  assume  equilibrium  equa¬ 
tions  of  state,  transport  processes,  and  reaction  rates,  are  not  valid. 
Typical  molecular  dynamics  simulations  (3)  describe  the  time  evolution  of 
the  dynamical  state  of  the  system  by  solving  the  classical  equations  of 
motion.  These  calculations  incorporate  all  of  the  interparticle  inter¬ 
actions  as  accurately  as  possible.  For  most  applications,  this  classical 
approximation  to  the  dynamics  is  justified  because  typical  quantum  mechani¬ 
cal  Interference  effects  (e . g. , resonant  cross  sections)  tend  to  be  averaged 
out  in  systems  possessing  many  degrees  of  freedom  (1). 

Manuscript  approved  October  1,  1985. 


'Jnder  some  circumstances  however,  there  might  be  residual  quantum 
mechanical  effects,  even  though  a  statistical  average  is  taken.  This  paper 
presents  the  preliminary  work  done  to  develop  phenomenological  models  of 
vibrational  energy  exchange  which  can  be  incorporated  into  large  scale 
classical  simulations.  These  models  mighc  be  designed,  for  examoLe,  to 
incorporate  some  of  the  effects  of  discrece  vibrational  excitations,  and 
thus  test  when  they  might  be  important  in  the  calculation.  The  first  step 
in  developing  such  models  is  a  careful  comparison  between  representative 
classical  and  quantum  calculations,  such  as  nresented  here. 

II.  Background  to  Present  Research 

Due  to  practical  computational  limitations,  most  molelcular  dynamics 
studies  treat  molecules  as  rigid  bodies  (5).  While  this  approach  precludes 
the  possibility  of  vibrational  energy  exchange,  it  does  incorporate  the  very 
important  dynamical  effects  of  rotational  and  translational  energy  exchange 
chough  molecular  collisions.  There  are  many  examples  in  the  literture  of 
simulations  involving  various  shapes  of  rigid  molecules  (6).  The  inter¬ 
action  between  molecules  is  generally  taken  as  a  sum  of  pairwise  additive, 
point  interactions  usually  (but  not  always  (7))  acting  between  atomic 
centers.  The  equations  of  motion  are  most  efficiently  integrated  in  a 
cartesian  coordinate  system  with  a  series  of  holonomic  constraints  imposed 
to  reoroduce  the  molecular  structure  (5). 

Some  dynamics  simulations  have  included  vibrational  degrees  of  freedom 
explicitly  (see  for  example,  reference  (9)).  This  requires  a  knowledge  of 
noc  onlv  the  intermolecular  interactions  (including  their  angular 
dependences),  but  also  a  knowledge  of  the  intramolecular  potential  as  a 
function  of  the  nuclear  coordinates.  Explicitly  including  the  vibrational 
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degrees  of  freedom  might  also  allow  for  the  additional  reactive  processes  of 
bond  dissociation  or  association  (10).  A  disadvantage  of  this  approach  is 
the  requirement  of  treating  all  degrees  of  freedom  on  the  same  dynamical 
footing,  t.e.,  classically.  While  the  classical  treatment  of  the  transla¬ 
tional  and  rotational  motions  is  generally  an  excellent  approximation,  the 
vibrational  interactions  begin  to  become  less  reliable  (8). 

A  further  disadvantage  of  explicitly  including  all  degrees  of  freedom 
in  a  molecular  dynamics  simulation  is  the  problem  of  temporally  resolving 
the  rapid  internal  vibrations.  This  requires  an  integration  time  step  on 
the  order  of  1/40  of  a  vibrational  period  for  systems  with  strong  but 
realistic  bond  potentials  (9).  Some  studies  have  avoided  this  problem  by 
including  artificially  soft  bonds,  which  vibrate  at  a  low  frequency  (10,11). 
The  less  rigorous  rigid-body  type  of  simulation  can  accomodate  a  larger  time 
step  in  the  numerical  solution  of  the  classical  equations  of  motion,  con¬ 
sequently  a  longer  period  of  time  can  be  simulated.  In  addition,  there  are 
generally  fewer  equations  of  motion  since  there  are  fewer  degrees  of  free¬ 
dom. 

The  goal  of  the  present  research  program  is  to  Investigate  ways  in 
which  the  effects  of  the  vibrational  degrees  of  freedom  might  be  included 
within  the  framework  of  the  computationally  tractable  rigid  body  approach  to 
molecular  dynamics  simulations.  Of  particular  interest  eventually  is  the 
determination  of  the  Influence  of  the  vibrational  degrees  of  freedom  on  the 
accumulation  of  internal  energy  and  subsequent  bond  dissociation.  The 
initial  step  in  such  a  study  must  be  a  determination  of  the  quality  of  a 
classical  or  semiclassical  treatment  of  a  vi brationaily  inelastic  molecular 
collision.  Such  a  determination  can  be  made  by  comparing  an  exact  quantum 
calculation  to  the  corresponding  classical  and  semiclassical  treatments. 
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Previous  work  in  comparing  quantum,  semiclassical  and  classical  methods 
for  creating  molecular  collision  phenomena  has  been  motivated  primarily  by 
interest  in  ascertaining  the  correspondence  between  classical  mechanics  and 
quantum  mechanics  (12).  The  motivation  for  the  present  work  stems  from  the 
need  to  incorporate,  within  the  classical  framework  of  dynamics  simulations, 
Che  important  effects  arising  from  vibrational  energy  exchange,  quantum 
effects  due  to  interference,  and  the  effects  arising  from  dissociation. 

Apart  from  the  computational  advantage  of  the  classical  approach  to 
molecular  dvnamics,  it  also  possesses  the  advantage  of  allowing  one  to 
follow  the  "trajectories,"  especially  through  the  important  interaction 
region.  Quantum  dynamical  treatments  allow  one  Co  "look"  at  the  physically 
measurable  quantities  only  before  and  after  collision,  preventing  intuitive 
insight  into  the  effect  of  a  third  body  or  other  perturbations  which  might 
interrupt  the  two  body  interaction.  Such  information  from  quantum  calcula¬ 
tions  would  be  extremely  useful,  and  attempts  have  been  made  previously  to 
address  this  question  (13). 

The  first  step  in  the  present  work  is  the  selection  of  the  particular 
model  molecular  collision  problem  to  serve  as  a  benchmark  for  investigation. 
The  model  chosen  involves  the  collinear  inelastic  collision  between  diatomic 
molecules,  the  details  of  which  are  described  in  section  III.  Transition 
probabilities  are  determined  quantum  mechanically,  semiclassically  (using  a 
classical  path  approximation)  and  classically.  The  precise  methodologies 
used  are  described  in  section  IV.  The  results  are  presented  and  compared  in 
section  7.  The  last  section  contains  a  discussion  of  the  results  and  how 
they  impact  on  the  question  of  improvement  of  molecular  dynamics  simula¬ 
tions. 
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III.  The  Model  Problem 


The  problem  we  consider  is  a  collinear  collision  between  two  homo- 
nuclear  diatomic  molecules.  Each  molecule  is  represented  as  a  harmonic 
oscillator.  The  masses  and  force  constants  are  chosen  to  correspond  to 
those  of  the  oxygen  and  nitrogen  molecules.  The  incermolecular  interaction 
arises  from  a  Lennard-Jones  ootential  acting  between  the  inside  atoms  in  the 
collinear  system.  Thus  only  adjacent  atoms  interact.  This  is  illustrated 
schematically  in  Fig.  I,  where  the  atoms  are  labeled  a,b,c  and  d  and  the 
potential  between  atoms  can  be  written: 

Vab(rl}  “  -r  (rl  -  rL0)2  (la) 

Vcd(C2)  ’  T  (r2  -  r20)2  (lb) 


((  y~  )12  -  (^2-  )6) 
'  be  be 


V 

ac 


0 


(Id) 


where  r^0  and  r2g  are  the  equilibrium  separations  of  the  molecules  ab  and  be 
respectively.  The  separation  between  the  centers  of  mass  of  the  two 
molecules  is  denoted  by  R, 


R  >  Rbc  m  1  /  2(  r  j_  +  r2  ) 


The  harmonic  force  constants  are  denoted  ’<  ^  and  ’<0 .  o  and  £  are  the  Leonard 
Jones  distance  ar  .  energy  oarameters  respectively.  For  the  problem 
considered  here, 
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r i 3  =  2.232  bohr 

r2  3  3  2.074  bohr 

ni  =  m,_  =  Lb  anw 
a  b 


'<3  =  20.57  av/(bohr)2 
k-,  *  40.09  ev/(bohr)2 
a  =•  4.730  bohr 
e  =■  .010  ev 


'•v'e  exnress  the  Lennard-Jones  ootential,  which  is  a  function  of  9,  ,  in  terras 

be 

of  the  be  seoaration  when  the  two  osciLlators  are  at  there  equilibrium 
positions  (F),  and  the  small  deviations  from  this  value  (x,  and  x2). 


F  =*  R  -  1/2  +  ^20^ 


(2a) 


:l  "  ri  “  r10 


(2b) 


x2  3  r2  ~  r20 

Using  Eons,  (2)  in  Eqn.  (lc)  the  Lennard-Jones  potential  becomes, 


Ubc(F,xt,x2)  -  4e  {(  77 T 


>12  -  (fTI 


F  -  j  (xL  +  x2) 


7  (xl  +  x2} 


)S> 


(3) 


4e  {(  ~  )12  (1  -  a)'12  -  (  £  )6  (1  -  a)'6} 


where, 


a  =  L/2(x,  +•  x2)/F  . 


(4) 


The  kinetic  energy,  in  the  center  of  mass  coordinate  system,  can  be 


written  as  , 


where  Oj,  o, ,  and  P  are  the  momenta  conjugate  to  r  ,  r7,  and  R  and  the 
reduced  masses  are  defined  as, 


ab  m^  +  m^ 


cd  m  +  m, 
c  d 


=  (TPa  *  V(mc  +  V 

ab,  cd  m  +  m.  +  m  +■  m , 
a  d  c  d 


The  Hamiltonian  for  Che  model  system  is  thus. 


P  P 

1  2 

(2u -  +  7  klxl2)  +  (2T“  +  7  k2x22> 


'2U  u  2 
ab 


+  W 

ab,cd 


ab,cd 


*  7U(!W 


IV.  Theoretical  Accroaches  to  the  Scattering  Problem 

A.  Quantum  Mechanical  Description 

Since  Che  Hamiltonian  for  the  system  does  not  depend  on  time,  we  search 
for  stationary  state  solutions  of  the  Schroedinger  equation.  In  the  quantum 
mechanical  treatment  of  inelastic  scattering,  it  is  convenient  to  work  with 
a  complete  set  of  functions  of  the  internal  coordinates  rt  and  r2.  The 
scattering  wavefunction,  which  is  a  function  of  Che  three  coordinates  r,  , 
and  R,  is  then  'written  as  a  linear  expansion  in  this  sec  of  functions,  where 
the  exoansion  coefficients  are  functions  of  R. 


■  (R,r- ,r2)  =  _  f  (R)*  (r,,r2) 
d  -  L  pn  n  1  ^ 
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The  sec  of  functions  (r  ,ri);  is  chosen  to  be  the  eigenfunctions  of  the 
internal  oortion  of  the  Harai Itonian.  Since  the  internal  portion  of  the 
Hamiltonian  is  the  sum  of  two  harmonic  oscillator  Hamiltonians,  the  eigen¬ 
functions  are  simply  products  of  harmonic  oscillator  eigenfunctions. 


H°  ( r  ,  r  -, )  i  (r  ,r7)  E  $  (r,  ,r,) 
n  4  -  n  n  4  - 


(9) 


whe  re 


&n^rl*r2)  =  ®iCr1)bj(r2) 


(10) 


and 


E  =*  E  +  Z. 
n  x  ] 


l  !<1  1/2  1  ,ft2  1/2 

(1  +  +  (j  + 

ab  cd 


(ID 


The  index  n  soecifies  the  state  of  both  molecules,  i.e.  specifies  i  and  j. 
The  time  independent  Schroedinger  equation  for  the  entire  system  is. 


r  -fO  d2 

'  2u  , 

ab,cd 


+  V, 


be 


(R,rI,r2)  H°  ( r  ]  ,  r  2 )  ]  l  f  (R) 

1  -  41  np 

n 


*  (r,. 
n  1 


a  E  ~  fon(R>  Vrl’r2>  <12) 

n 

hultiolying  on  the  left  bv  ,r,)  and  integrating  over  the  internal 

coordinates,  r,  and  r7  yields  the  following  set  of  coupled  equations  for  cl 
exDansion  coefficients  f  ^(R), 


2u 


ab,cd 


,00 


(R)f  (R) 


03) 


where 


VR) 


fti(rl*r2)Vbc(R,ri,r2)  *a(rl>r2)dridr2 


(14) 


Eqn.  (13)  can  be  rewritten  in  matrix  form  as, 


dR 


•2  £(R)  =  D(R)f(R) 


(L5) 


with  the  symmetric  matrix  D  defined  by, 

2U 


Du  •  Djt  •  'VR)  *  (v  E>V 


(16) 


The  coupled  equations  (Eq.  13)  are  solved  numerically  by  the  R-matrix  method 
(14).  The  essential  idea  of  the  method  is  to  solve  the  problem  in  the 
interaction  region  where  the  basis  states  are  coupled  to  one  another  and 
then  to  match  this  wavefunction  with  the  known  uncoupled  asymptotic  forms. 
The  solution  in  the  interaction  region  is  determined  by  dividing  the  R-space 
into  many  sectors  within  which  the  coupling  potential  does  not  change 
significantly.  Within  each  sector,  the  unitary  transformation  is  found 
which  diagonalizes  the  coupling  matrix  D. 


U+(R)  D(R)U(R)  =  D'(R) 


(17) 


Equation  (15)  may  then  be  written  within  each  sector  in  the  uncoupled  form, 


where 


d2  £'(R)  -  D'f(R)  (18) 

dR2 
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f  ' (R)  =  0  (R)£(R) 


(19) 


and  D'  is  a  diagonal  matrix.  The  solacion  is  Chen  snatched  ac  adjacent 
sectors,  ensuring  the  continuity  of  the  wavefunction  and  its  derivative. 

Formallv,  the  summations  over  n  in  Eons.  (12)  and  (13)  extend  to 
infinity.  This  is  because  an  •'nfinita  number  of  the  functions  5n(r,  »r2)  are 
required  to  form  a  complete  set  of  functions  which  soan  the  space  of  r.  and 
r2 .  The  choice  of  internal  state  eigenfunctions  to  represent  this  space  is 
a  judicious  one  however,  since  it  is  expected  that  even  in  the  interaction 
region,  the  oscillators  will  to  some  extent  retain  their  identity.  In 
practice  we  do  a  finite  basis  set  expansion.  Functions  are  selected  on  the 
basis  of  their  energy,  and  the  calculations  are  performed  until  increasing 
the  length  of  the  expansion  has  no  effect  on  the  answer. 

To  solve  the  Schroedinger  Equation,  we  need  to  know  the  elements  of  the 
coupling  matrix  defined  in  Eqn.  (14)  as  a  function  of  the  molecular  separa¬ 
tion  R.  This  requires  either  the  calculation  of  N-squared  integrals  (N  is 
the  number  of  basis  states)  each  time  a  different  value  of  R  is  considered 
during  the  course  of  the  calculation,  or  the  prior  storage  of  a  large  amount 
of  R-dependent  data.  To  avoid  the  storage  problem  and  to  simplify  the  inte¬ 
gral  computation,  we  use  the  following  scheme  for  evaluating  the  matrix 
elements  of  the  Leonard- Jones  potential. 

In  Eqn.  (3),  the  Lennard-Jones  potential  acting  between  the  inside  atoms 
was  written  so  as  to  separate  the  dependence  of  the  Dotential  on  the  inter¬ 
nal  coordinates  (the  a  dependence)  from  the  dependence  on  the  collision 
coordinate  (the  F  dependence).  For  a  given  molecular  seoaration,  we  can  now 
exoand  the  potential  about  the  equilibrium  positions  of  the  oscillators. 

Thus  equation  (3)  is  expanded  in  a  Taylor  series  in  a  about  the  point  n=*0. 


10 


(20) 


Carrying  the  expansion  to  fourth  order  and  integrating  over  the  internal 
coordinates,  the  coupling  matrix  V  becomes, 


The  important  point  is  chat  these  matrices,  and  hence  the  integrals  which 
make  up  the  coupling  elements,  have  no  dependence  on  the  collision  coordi¬ 
nate.  The  coupling  matrix  is  obtained  as  a  function  of  R  by  multiplying 
these  previously  stored  matrices  by  the  R-dependent  coefficients  shown  in 
Sqm.  (21).  We  have  found  that  Che  coupling  matrix  is  well  converged  when 
the  expansion  of  the  Lannard-Jones  potential  is  carried  out  to  fourth  order. 
As  an  additional  test,  we  have  performed  a  number  of  classical  trajectory 
calculations  using  thi3  truncated  expansion  of  the  potential.  Again,  the 
fourth  order  expansion  produces  results  which  are  the  same  as  the  analytical 
form  for  the  potential  to  four  decimal  places. 
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3.  Semiclassical  Description 

The  semiclassical  method  we  have  chosen  is  generally  referred  to  as  a 
classical  path  approach.  The  underlying  idea  is  that  motion  along  a 
collision  coordinate  is  responsible  for  initiating  transitions  among  the 

discrete  quantum  states  of  the  colliding  partners,  but  that  this  motion  is 
essentially  classical  in  nature.  One  first  chooses  a  reasonable  mean  tra¬ 
jectory  for  the  collisional  degree  of  freedom.  This  trajectory  could  be,  for 
example,  a  classical  trajectory  associated  with  a  hypothetical  system  for 
which  the  internal  degrees  of  freedom  have  been  frozen.  Since  the  inter¬ 
action  among  the  collision  partners  is  a  function  of  the  collision  coordi¬ 
nate,  this  predetermined  trajectory  provides  the  interaction  potential  as  a 
function  of  time.  Using  this  interaction  potential,  one  then  constructs  a 
time  dependent  Hamiltonian  operator  for  the  internal  degrees  of  freedom. 

The  time  dependent  Schroedinger  equation  is  then  solved  in  the  subspace  of 
the  internal  degrees  of  freedom. 

For  our  model  problem,  we  first  separate  the  Hamiltonian  into  a  portion 
which  contains  all  of  the  internal  coordinate  dependence  and  a  portion  which 
depends  upon  the  collision  coordinate  only, 

2  2 

H  -  +  V(R,r  )}  +  {H°  +  V(R,r)  -  V(R,r  )}  (23) 

ZU  ,  ,  dR.  0  o 

ab,cd 

The  first  set  of  brackets  contain  a  simple  Hamiltonian  for  the  collision 
coordinate.  We  solve  this  one  dimensional  problem  classically  to  obtain  the 
collision  coordinate  R  as  a  function  of  time.  The  second  set  of  brackets 
contain  a  Hamiltonian  for  the  internal  degrees  of  freedom.  The  collision 
coordinate  appears  as  a  oarameter  in  this  Hamiltonian,  its  time  dependence 
having  been  pre-determined.  The  equation  which  must  be  solved  is  then 
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3'?(r  L  ,r2,c) 


[H°  +•  V(R,r)  -  V(R,r  )}  »(r,,r,,t)  (24) 

0  c 


As  in  che  close  counting  calculations,  we  expand  Che  wavefunction  in  a  basis 
of  eigenstaces  of  H°  ,  che  harmonic  oscillator  produce  functions. 


Left  multiplying  Eqn.  (24)  by  one  particular  eigenfunction  and  integrating 
over  the  internal  coordinates  produces  equations  describing  the  time  depen¬ 
dence  of  che  real  and  imaginary  components  of  the  expansion  coefficients: 


(E.  -  V(R,r  )}s.  + 
J  03 


t!j  -  - 1  v* 


(26) 


These  equations  can  then  be  solved  numerically.  A  typical  calculation 
begins  with  the  a  large  separation  between  the  two  molecules  with  unit 
probability  of  being  in  one  of  the  internal  eigenstates  (entrance  channels). 
The  system  then  proceeds  through  the  interaction  region  and  back,  out  to 
large  separation.  The  square  moduli  of  the  expansion  coefficients  then  give 
the  probability  that  a  measurement  will  find  the  system  in  each  of  the 
possible  exit  channels.  The  validity  of  the  classical  path  approach  is 
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discussed  by  Child  [15].  Essentially  what  is  found  is  what  one  intuitively 
expects.  The  change  in  momentum  associated  with  transitions  among  various 
diabatic  states  should  be  small  relative  to  the  momentum  associated  with  the 
collision  coordinate.  In  other  words,  the  classical  paths  associated  with 
different  entrance  and  exit  channels  should  all  lie  close  to  the  assumed 
oath. 


C.  Classical  Description 

The  classical  descriotion  of  the  inelastic  scattering  is  also  based  on 
the  Hamiltonian  of  Eqn.  (7).  While  classical  trajectories  governed  by  this 
Hamiltonian  are  straightforward  to  calculate  (16),  it  is  less  clear  how  to 
compare  the  trajectory  results  with  the  corresponding  quantum  mechanical 
results.  The  most  widely  applied  approach  for  making  comparisons  between 
classical  and  quantum  results  is  the  quasi-classical  method  (17).  Before 
applying  this  approach,  a  canonical  transformation  (26)  is  applied  to  the 
Hamiltonian  of  Eqn.  (7)  such  that  the  new  internal  coordinates  and  momenta 
correspond  to  the  action-angle  variables.  The  action  variables  are 
proportional  to  the  total  energy  of  the  oscillators  and  are  therefore 
susceptible  to  being  "quantized"  in  units  of  Planck's  constant  (18).  The 
angle  variables  then  correspond  to  the  phase  of  the  oscillators.  For  the 
harmonic  oscillators  under  consideration,  the  old  internal  coordinates  and 
momenta  (r,p)  are  given  in  terms  of  the  new  action-angle  variables  (n,o)  as: 


r  -  r 

o 


r  2( n+  1/2)1*, 1/2  , 

[ - - - ]  sin  q 


p  =  [2(n  +  l/2)fiuw]i/2  cos  q 


(27) 
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The  Hamiltonian  of  Eqn.  (7),  in  terms  of  these  new  conjugate  variables,  is 
given  by, 


H(P,R,n,n2  ,qi,q2)  =*  -y- -  +V  (R,r1  Cn  ,  ,  q  1 ) ,  r2  (n2  ,  q2 ) ) 

ab.cd 


*  (n  i  *  i/2)  fi'H.  +  (n,  *  1/2) 


A  useful  property  of  the  internal  action-angle  variables  is  that  in  the 
asymototic  region,  where  the  interaction  potential  vanishes,  the  Hamiltonian 
(Eqn.  (28))  becomes  independent  of  the  new  angle  variables.  Hamilton's 
equations  of  motion,  given  by. 


7n\  -  ql 


3H 

lq,  *  nl 


JKZ  M  q2 


therefore  imply  that  the  action  variables  (i.e.,  the  internal  oscillator 
energies)  become  constants  of  the  motion  both  before  and  after  collision. 

The  major  obstacle  in  making  comparisons  between  classical  and  quantum 
mechanical  results  is  the  fact  chat  the  classical  action  variables  (i.e., 
the  classical  counteroart  of  the  vibrational  quantum  numbers)  are  allowed  to 
have  a  continuous  range  of  values,  whereas  the  quantum  mechanical  oscilla¬ 
tors  are  allowed  to  have  only  discrete  values  for  quantum  numbers  (i.e., 
integers).  The  crux  of  the  quasi-classical  aooroach  (also  known  as  the 


y\'> 


•  ’  /  ■vAv/>  i 
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"binning”  or  "histogram"  method)  is  simply  to  assign  final  (asymptotic)  non- 
integer  values  of  the  action  variables  to  the  nearest  integer  value.  In 
practice,  a  large  number  of  trajectories  are  calculated,  each  trajectory  is 
subsequently  "binned",  and  a  statistical  analysis  is  performed  on  the 
results  of  binning  the  "batch"  of  trajectories. 

Trajectories  begin  in  the  asymototic  region,  where  the  variables  are 
assigned  the  following  initial  conditions  (assuming  a  total  collision  energy 
5): 

R  *  large 

nt  *  integer  (corresponding  to  initial  quantum  state) 

a,  *  integer  (corresponding  to  initial  quantum  state)  (30) 

P  *  [21Jab,cd(E  "  (nl  +  *  C«2  +  1/ 2>Ti^2)  ] 1  /2 

q  j_  ■  random  (0,2m) 

q2  *  random  (0,2m). 

The  angle  variables  (i.e.,  the  initial  oscillator  phases)  are  chosen 
randomly  from  a  uniform  distribution  between  0  and  2m.  A  large  number  of 
trajectories  simulates,  to  a  close  approximation,  all  of  the  possible  types 
of  molecular  encounters  for  this  collinear  model.  The  equations  of  motion 
(Eqn.  (29))  are  numerically  integrated  until  the  collision  coordinate 
becomes  large.  At  this  point,  the  final  action  variable  is  assigned  to  a 
"bin"  corresponding  to  the  quantum  final  number. 
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V.  Computational  Results 

Calculations  were  performed  for  total  collision  energies  ranging  from 
1.0  ev  to  2.25  ev.  Figure  II  shows  the  intermolecular  potential  acting 
between  atoms  b  and  c.  On  the  same  energy  scale,  the  figure  also  shows  the 
energy  level  soacings  for  the  two  oscillators  and  the  oositions  of  the  first 
several  levels  of  the  combined  system,  'fate,  for  example,  that  a  collision 
starting  in  the  ground  state  (0,0)  with  a  total  energy  of  1.0  ev  would  have 
roughly  .21  ev  as  initial  internal  (vibrational)  energy  and  would  thus  have 
about  .76  ev  as  initial  translational  energy  in  the  collision  coordinate. 

The  Lennard-Jones  distance  parameter,  a  =*  2. 5  A,  is  somewhat  smaller  than 
the  approximately  3.5  A  that  one  might  expect  to  see  for  a  nonbonded  inter¬ 
action  between  nitrogen  and  oxygen.  This  was  done  to  increase  the  amount  of 
inelastic  scattering  at  low  energies.  The  indexing  scheme  for  the  diatomic 
product  functions  is  shown  in  Table  I. 

As  was  mentioned  earlier,  the  flexibility  in  the  semiclassical  approach 
is  in  the  choice  of  the  classical  path.  A  strict  interpretation  of  a 
quantum  mechanical  description  of  the  scattering  event  does  not  recognize 
the  existence  of  well  defined  "path"  for  the  collision  coordinate  along 
which  the  system  proceeds.  We  do,  however,  speak  of  such  a  path  as  it 
relates  to  the  approximate  semiclassical  description.  As  described  above, 
the  most  obvious  choice  for  this  path  is  a  classical  trajectory  for  the 
collisional  degree  of  freedom  with  the  internal  coordinates  frozen  at  their 
equilibrium  values. 

Table  IT  displays  a  comparison,  for  a  particular  collision,  of  quantum 
mechanical  and  semiclassical  state-to-state  probabilities  where  various 
classical  oath  approaches  were  used.  The  semiclassical  calculations  at  this 
energy  all  use  25  basis  states.  The  close  coupling  calculations  use  30.  At 
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this  energy,  14  of  these  channels  are  ooen.  The  total  energy  for  tne  calcu¬ 
lation  is  1.25  ev.  and  the  oscillators  are  initially  in  their  ground  states. 
At  1.25  ev,  all  of  the  calculations  show  roughly  90  percent  retention  prob¬ 
ability  (elastic  scattering). 

The  first  semiclassical  column  in  Table  II  corresponds  to  a  calculation 
in  which  the  oath  is  the  simple  classical  trajectory  for  the  collision 
coordinate.  The  agreement  is  at  least  qualitatively  quite  good.  'tote  the 
direction  of  the  deviation.  The  semiclassical  scheme  predicts  too  much 
inelastic  scattering.  This  deviation  can  be  understood  by  noting  that 
energy  is  not  conserved  in  this  calculation.  The  system  begins  with  the 
minimum  allowable  energy  in  the  internal  vibrations.  As  the  molecules 
proceed  along  the  classical  trajectory,  the  Hamiltonian  for  the  internal 
degrees  of  freedom  changes  with  time  reflecting  the  coupling  between  the 
oscillators.  This  coupling  produces  an  amplitude  to  find  the  system  in 
higher  energy  internal  states.  Energy  is  not  a  constant  of  the  motion  for 
the  internal  portion  of  the  motion  since  the  Hamiltonian  depends  on  time, 
but  energy  is  conserved  by  construction  for  the  collision  coordinate. 
Therefore  energy  is  not  conserved  for  the  system  as  a  whole.  The  molecules 
can  emerge  from  the  collision  with  appreciable  probabilities  of  being  found 
in  higher  energy  states,  yet  the  molecules  translate  away  from  one  another 
with  the  same  kinetic  energy  that  they  had  during  the  approach.  There  is  a 
net  edition  of  energy  to  the  system. 

In  the  collision  considered  here,  a  consequence  of  this  energy  addition 
is  too  much  inelastic  coupling.  .As  the  molecules  approach  each  ocher  and 
translational  energy  begins  to  be  transferred  to  internal  vibracional 
energv,  one  expects  the  molecules  to  slow  down.  With  less  energy  in  the 
translation,  the  molecules  would  not  approach  as  closely  and  consequently, 


18 


'here  would  be  less  coupling.  A  prediction  of  the  probability  for  a  system 
starting  in  a  highly  excited  state  to  emerge  in  the  ground  vibrational  state 
has  the  opposite  problem.  Starting  with  only  a  small  portion  of  the  energy 
in  translation,  the  trajectory  would  not  provide  enough  coupling  since  it 
cannot  reflect  a  transfer  of  energy  from  the  internal  coordinates  to  the 
collision  coordinate. 

These  two  probabilities,  of  a  transition  from  the  ground  state  to  a 
particular  excited  state  and  from  that  excited  state  to  the  ground  state, 
are  known  to  be  equal  by  microscopic  reversibility.  It  is  then  reasonable 
to  symmetrize  this  probability  matrix  to  impose  this  feature  of  the  correct 
ohysical  result.  We  do  this  simply  by  taking  the  arithmetic  mean  of  two 
symmetrically  disposed  probabilities.  The  diagonal  elements,  the  probabili¬ 
ties  of  elastic  scattering,  are  Chen  adjusted  by  the  normalization  condi¬ 
tion.  Results,  which  have  been  modified  in  this  way  are  shown  in  column  b 
of  Table  II.  These  are  a  notable  improvement  (in  comparison  to  the  quantum 
mechanical  calculations)  over  the  probabilities  in  the  first  column. 

The  problem  of  the  lack  of  energy  conservation  in  semiclassical  calcula¬ 
tions  of  this  sort  is  well  known.  Schemes  designed  to  compensate  for  this 
weakness  have  been  proposed,  usually  by  choosing  a  classical  trajectory 
which  is  characterized  by  an  energy  (19)  or  velocity  (20)  averaged  between 
initial  and  final  states.  These  schemes  yield  a  path  which  is  suitable  for 
the  one  state-to— state  transition  of  interest.  We  are  interested  in  simul¬ 
taneously  ascertaining  a  orobability  of  the  system  emerging  in  each  of  the 
exit  channels  given  that  it  began  in  a  particular  entrance  channel.  The 
above  discussion  suggests  the  incorporation  of  some  sort  of  feedback  mechan¬ 
ism  in  the  semiclassical  calculations.  Perhaps  the  classical  oath  can  be 
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dynamically  modified  so  as  Co  reflect,  at  least  partially,  the  vibrational- 
translational  energy  transfer.  Ideas  along  these  lines  have  been  prooosed 
in  the  past,  focusing  on  the  coupling  of  the  oscillator  response  to  the 
classical  path  (21)  or  on  the  incorooration  of  an  effective  potential  for 
the  collision  coordinate  (22). 

The  probabilities  in  the  third  semiclassical  column  of  table  II,  column 
c,  result  from  a  simnle  incorporation  of  some  "back  coupling"  of  the  inter¬ 
nal  degrees  of  freedom  to  the  relative  motion  of  the  collision  partners. 

The  classical  path  at  the  onset  of  the  collision  is  the  unperturbed  classi¬ 
cal  trajectory  associated  with  the  total  collision  energy  less  the  internal 
energy  for  the  particular  entrance  channel.  The  path  is  modified  at  each 
timestep  of  the  numerical  intregration  of  San.  (26)  according  to  an  energy 
conservation  constraint  imposed  on  the  system  as  a  whole.  Since  the  time 
sceo  can  be  made  arbitrarily  small,  this  corresponds  to  a  continuous  modifi 
cation.  The  total  energy  at  any  instant  is, 


E  *  constant 


2u 


Y  Y  E  +  V 

*i  h  i  int 


ab,cd 


(31) 


The  first  and  second  terms  are  the  kinetic  energy  in  the  collision  coordi¬ 
nate  and  the  internal  energy  respectively.  The  third  term  is  the  inter¬ 


action  potential  which  is  calculated  as, 


The  kinetic  energy  in  the  collision  coordinate,  and  hence  the  velocity  are 
continuously  modified  such  that  Fqn.  (31)  remains  satisfied.  The  relative 
motion  is  essentially  governed  by  an  effective  potential  which  is  a  function 
of  the  Quantum  mechanical  state  of  the  internal  system  as  well  as  the  rela¬ 
tive  center  of  mass  seoaration.  The  resulting  path,  although  not  ootimal 
for  any  oarticular  state-to-state  probability  is  biased  toward  exit  channels 
which  olav  the  greatest  role  in  the  inelastic  scattering. 

The  results  in  column  c  of  Table  II  show  that  modifying  the  classical 
oath  as  described  above  improves  the  orobabilities  compared  to  case  a, 
almost  halving  the  errors  in  this  case.  Symmetrizing  the  probability  matrix 
as  before  leads  to  the  entries  found  in  the  last  column  of  Table  II.  These 
results  are  in  excellent  agreement  with  the  quantum  mechanical  probabili¬ 
ties.  All  of  Che  semiclassical  results  reported  below  have  been  obtained  by 
this  procedure  of  imposing  energy  conservation  and  symmetrizing.  Although 
the  agreement  between  the  quantum  mechanical  and  semiclassical  results  is 
not  always  quite  as  striking  as  it  is  in  the  above  example,  we  find  that  it 
remains  very  good  for  every  case  we  studied.  Figure  III  shows  inelastic 
collision  probabilities  for  the  conditions  in  the  sample  calculation  dis¬ 
cussed  above.  Shown  are  probabilities  that  a  system  starting  in  a  particu¬ 
lar  initial  state,  has  a  collision  which  is  inelastic,  leading  collectively 
to  any  ocher  state.  That  is,  we  are  plotting  one  minus  the  retention  prob¬ 
ability.  The  close  couoled  quantum  results  are  labeled  with  a  "Q"  and  the 
semiclassical  results  are  labeled  with  an  "S”.  The  classical  results,  also 
shown  and  labeled  with  a  "C,"  ’will  be  discussed  later.  Although  there  are 
fourteen  open  channels  at  this  energy,  the  probability  of  an  inelastic  event 
approaches  essentially  zero  by  about  the  ninth  or  tenth  state.  Figure  IV 
corresponds  to  collisions  at  1.0  ev,  a  lower  total  energy  than  Figure  III. 
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Figure  V  corresponds  to  a  Higher  ener^  1.75  ev.  Mote  "hat  Che  scale  of  Che 
graph  has  been  expanded  by  an  order  of  magnitude  for  Che  1.0  ev  collisions 
due  Co  Che  small  orobabilicy  of  inelascic  scaccering.  For  the  1.0  ev  calcu¬ 
lations,  there  are  eight  open  channels.  Sixteen  basis  states  were  used  in 
the  quantum  calculations  and  fifteen  in  the  semiclassical  calculations.  As 
seen  in  Figure  V,  there  is  considerably  more  inelastic  scattering  at  1.75 
ev,  where  forty  five  basis  states  were  used  for  both  the  quantum  and  the 
semiclassical  calculations.  Figures  III-V  show  basically  the  same  trends  in 
the  probabilities.  Increasing  the  energy  simply  increases  the  amount  of  the 
inelascic  scattering. 

It  is  also  instructive  to  compare  the  quantum  and  semiclassical  results 
for  scaccering  from  a  particular  initial  state.  Figure  VI  shows  such 
results  for  a  collision  at  1.75  ev  starting  in  the  second  state,  i.e.,  a 
state  with  Che  nitrogen  molecule  in  its  ground  state  and  the  oxygen  molecule 
in  its  first  excited  state.  .As  shown  in  Figure  V,  this  is  the  initial  state 
for  which  the  semiclassical  and  quantum  calculations  show  the  least  agree¬ 
ment.  There  is  roughly  equal  probability  for  a  transition  to  the  first 
state  and  a  transition  to  the  fourth  state.  These  correspond  to  transitions 
of  the  oxygen  to  its  ground  state  and  to  its  second  excited  state  respec¬ 
tively.  The  semiclassical  calculations  slightly  underestimate  this  prob¬ 
ability  in  both  cases.  The  normalization  condition  requires  chat  the  sum  of 
the  deviations  in  Figure  VI  vanish.  It  therefore  follows  that  Che  retention 
probability  is  overestimated. 

The  more  complete  results  of  Che  semiclassical  and  quantum  mechanical 
calculations  at  the  three  energies  (1.0  ev,  1.25  ev,  and  1.75  ev)  are  shown 
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In  Tables  III-V.  The  serai classical  stace-to-state  probabilities  are  dis- 
olaved  in  parentheses  below  the  quantum  probabilities.  The  quantities  plot¬ 
ted  in  Figures  III,  IV,  and  V  are  simply  one  minus  the  diagonal  elements  in 
these  Tables.  Figure  VI  is  the  second  column  in  Table  V.  Not  all  of  the 
twentv  seven  available  channels  are  represented  in  Table  V.  One  can  see 
however  that  che  amount  or  inelastic  scattering  is  dropping  off  rapidly. 

Tie  classical  results  shown  in  figures  III-K  are  all  based  on  the 
auas i-class ical  (17)  (binning)  method  described  in  section  III.  Transition 
probabilities  from  each  specified  initial  state  were  calculated  from  500 
trajectories.  Previous  applications  of  the  quasi-classical  method  have  had 
greatest  success  when  many  final  states  are  dynamically  accessible  (23).  It 
seems  reasonable  to  expect,  therefore,  that  higher  collision  energies  should 
begin  to  show  better  agreement  between  the  classical  results  and  the  corre¬ 
sponding  quantum  results.  3y  inspection  of  Figures  III-V  and  also  from  the 
results  at  2.25  ev  shown  in  Figure  VII,  one  can  see  that  the  agreement  does 
not  improve  as  expected.  Certain  state-to-state  probabilities  do  show 
excellent  agreement,  however,  and  general  trends  in  transition  probabilities 
are  modelled  reasonably  well.  For  example,  for  collisions  at  1.75  ev.  (Fig. 
VI),  the  classical  estimate  of  the  elastic  scattering  from  the  n(02)  a  1, 
n(N2)  ■  0  initial  state  is  significantly  less  than  the  quantum  prediction, 
giving  rise  to  a  classical  overestimate  of  the  inelastic  scattering  from 
that  state  (Fig.  V),  At  2.25  ev.  (Fig.  VIII),  the  elastic  scattering  from 
the  ground  vibrational  state  is  in  excellent  agreement  with  the  quantum 
results.  At  2.25  ev  (the  highest  energy  considered),  the  classical  estimate 
is  high  for  inelastic  scattering  from  initial  staces  where  nitrogen  is  more 
excited  chan  oxygen,  and  it  is  low  for  initial  states  where  oxygen  is  more 
excited  than  nitrogen.  This  trend  Is  not  apparent  at  lower  collision 
energies . 
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At  low  collision  energy  (1.00  ev.),  the  classical  binning  aporoach 
shows  essentially  zero  inelastic  scattering  for  all  initial  states.  This  is 
a  severe  test  of  the  binning  procedure  and  it  is  normally  under  these  low 
energy  conditions  that  other  classical  aporoaches  e.g.,  the  moment  method 
(24)  or  Wigner  distribution  methods  (25)  have  had  more  success.  The 
unbinned  final  actions  for  a  set  of  400  trajectories  are  shown  in  Figure 
7TII.  It  is  clear  that  the  artificial  imoosition  of  binning  "boundaries" 
may  contribute  to  the  inadequacy  of  the  classical  results.  Other  procedures 
(such  as  those  mentioned)  may  prove  more  useful.  Work  is  in  progress  to 
determine  the  success  and  practicality  of  these  other  approaches. 

Similar  final  action  plots  are  shown  in  Figures  IX  and  X  for  a  higher 
collision  energy  (2.25  ev),  where  inelastic  scattering  dominates.  The 
initial  classical  actions  for  the  collisions  represented  in  Fig.  IX  corres¬ 
pond  to  quantum  state  8  and  for  Fig.  X  to  quantum  state  13.  The  energetic¬ 
ally  accessible  region  for  this  high  energy  system  is  extensive.  There  are 
44  open  channels  at  this  energy.  This  means  that  at  least  44  of  the  bins 
represented  in  Figures  IX  and  X  are  accessible  within  the  constraint  of 
energy  conservation.  The  emerging  boundaries  of  the  points  in  Figures  IX 
and  X  show  that  the  dynamically  accessible  regions  are  much  more  restricted. 
Figure  IX  for  example,  only  14  channels  are  being  populuated  (21  channels 
are  oooulated  via  the  corresponding  quantum  mechanical  dynamics) .  Apparent¬ 
ly,  the  classical  dynamics  precludes  the  possibility  of  populating  some  of 
these  bins  which  correspond  to  states  which  are  quantum  mechanically  access¬ 
ible.  Perhaos  an  alternative  choice  of  initial  condition  selection  (e.g., 
via  the  Wigner  distribution  method  (25))  or  analysis  of  final  conditions 
(e.g.,  the  moment  method  (24))  would  eliminate  some  of  the  discrepancies . 


VI.  Summary  and  liscussion 

As  a  shock  passes  through  a  condensed  phase  material,  each  molecule 
instantaneously  feels  an  impulse.  This  impulse  is  the  sum  of  the  changes  in 
the  individual  forces  the  molecule  feels  as  each  of  its  neighbors  moves. 
Studying  the  microscopic  behavior  of  such  a  system  requires  evaluating  how  a 
soecific  impulse  effects  a  molecule,  given  the  molecule's  initial  state  and 
the  impulse  characteristics.  To  address  this  oroblem,  we  first  need  to  ask 
about  the  behavior  of  the  microscopic  system,  and  then  about  a  macroscopic 
ensemble  of  many  such  systems. 

This  paper  addressed  part  of  the  miscroscopic  question.  Specifically, 
we  asked  how  quantum  mechanical  effects  influence  the  description  of  a 
microscopic  collision,  and  what  is  a  good  calculation  of  the  collision 
properties.  We  considered  two  harmonic  oscillators  interacting  through  a 
Lennard-Jones  potential.  This  model  problem  is  a  vehicle  for  studying  the 
quantum  effects  of  discrete  vibrational  states  on  the  collision  of  two 
molecules.  The  model  problem  was  solved  quantum  mechanically  by  a  close¬ 
coupling  method,  semiclassically  by  several  variations  of  a  classical  path 
approach,  and  classically  by  a  quasiclassical  trajectory  method. 

The  quantum  mechanical  calculation  was  considered  the  correct  answer 
against  which  we  compared  the  results  of  the  classical  and  semiclassical 
calculation.  The  output  of  this  calculation  was  the  final  distribution  of 
internal  staces  of  the  two  molecules,  once  they  had  collided  and  completely 
separated.  The  quantum  mechanical  results  were  discussed  in  Section  V. 

The  straightforward  semiclassical  classical  path  description  discussed 
in  Section  V  gave  reasonable  answers.  However,  when  the  path  was  modified 
to  conserve  energy  and  the  microscopic  reversibility  was  imposed  on  the 
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final  scace  Co  state  probabilities,  the  results  were  found  to  be  in 
excellent  agreement  with  the  close-coupling  calculations.  The  method 
faithfully  reproduced  the  exact  quantum  mechanical  stace-to-seate  transition 
probabilities  for  a  wide  range  of  collision  energies.  This  range  extends 
from  energies  which  produce  almost  no  inelastic  scattering  to  energies  which 
produce  mostly  inelastic  scattering.  This  result  is  encouraging  since  an 
attempt  was  made  to  choose  potentials  which  correspond  to  what  might  he 
expected  in  a  realistic  molecular  collision. 

The  classical  method  used  in  this  paper  is  the  quasi -classical 
trajectory  method.  Final  values  of  internal  coordinates  (action  variables) 
were  assigned  to  quantum  states  by  a  simple  binning  procedure.  The  results 
for  the  final  distribution  of  vibrational  states  of  the  molecules  do  not 
agree  particularly  well  with  the  quantum  and  semi -classical  calculations. 

In  particular,  the  region  of  phase  space  required  in  order  to  populate 
states  which  are  populated  quantum  mechanically  and  semiclassically  appears 
in  some  cases  to  be  dynamically  inaccessable. 

It  is  reasonable  to  ask  whether  modifying  the  classical  approach  would 
give  better  answers.  For  example,  the  quasi-classical  method  as  implemented 
is  not  microscopically  reversible.  This  important  dynamical  concept  may  be 
incorporated  in  the  calculations  by  selecting  the  initial  action  variables 
from  a  uniform  distribution  centered  around  a  particular  quantum  number. 

Such  a  microscopically  reversible  quasi-classical  technique  was  used  in  a 
few  of  the  initial  states,  but  it  showed  no  significant  effect  on  the 
results.  This  is  net  surprising  because  the  original  batches  of 
trajectories  lid  not  violate  microscopic  reversibility  to  any  significant 
degree.  Cther  approaches  to  final  action  analysis,  such  as  the  moment 
method,  or  a  modified  choice  of  initial  conditions  which  incorporate  quantum 


effects  more  directly,  such  as  the  Vigner  distribution  method,  might  improve 


the  agreement  between  classical  and  quantum  results. 

The  encouraging  result  of  this  paper  is  the  good  agreement  between  the 
quantum  and  the  semiclassical  predictions.  This  introduces  an  interesting 
possibility  for  molecular  dynamics  calculations.  The  semiclassical  calcula¬ 
tions  could  be  used  to  determine  the  distribution  of  states  arising  from  a 
collision.  Then  this  distribution  could  be  incoroorated  as  a  submodel  in  a 
classical  molecular  dvnamics  calculation.  This,  however,  involves  several 
leaps  of  faith.  First,  there  is  a  basic  problem  that  the  quantum  calcula¬ 
tion  only  gives  the  final  distribution  of  states,  in  the  asymptotic  region 
of  the  collision.  It  is  not  at  all  clear  how  well  the  semiclassical  calcu¬ 
lations  represent  the  collision  in  the  interaction  region.  The  premise  of 
the  molecular  dynamics  calculation  would  be  that  the  results  of  an  impulse 
felt  by  a  molecule,  due  to  its  simultaneously  interacting  with  many  mole¬ 
cules,  could  be  related  to  the  results  of  an  impulse  of  the  same  magnitude 
felt  by  the  collision  with  just  one  molecule.  Even  if  the  answers  are  not 
equivalent,  does  statistical  averaging  make  them  better.  These  are  just  a 
few  of  the  questions  that  will  be  addressed  in  the  future. 
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TABLE  II 


Comparison  of  quantum  mechanical  state  to  state  transition  probabilities 
with  those  calculated  using  a  variety  of  semiclassical  classical 
path  schemes  for  a  collision  with  total  energy  1.25  ev.  starting 
in  the  ground  state. 


STATE 


QUANTUM 


SEMICLASSICAL 


TABLE  III 


Quantum  mechanical  state  to  state  transition  probabilities  for 
model  system  with  a  total  collision  energy  of  1.00  ev.  Semiclassical; 
probabilities  (see  text)  are  shown  in  parentheses.  There  are  eight 


open  (energetically  accessible) 
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(.022) 
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TABLE  IV 


M 

[I 


Same  as  TABLE  III  with  a  total  collision  energy  of  1.25  ev.  There 
are  fourteen  open  channels. 
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( 

000) 

(.300) 

(.000) 

(.004) 

(.000) 

( 

300) 

(.995) 

(.000) 

( 

001) 

( 
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12 
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.000 
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.000 
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.300 
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( 
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( 
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( 
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13 
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.000 
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.000 

.000 
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.000 

.000 
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1 
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.300 
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( 
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( 

000) 
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( 
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(.000) 

( 
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(1 
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(.000) 

14 

.300 

.000 

000 
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.000 

.000 

.000 

.  000 

300 

.000 

.000 

000 

300 

1.00 

(.000) 

(.000) 

( 

300) 

( 

000) 

(.000) 

(.000) 

(.000) 

(.000) 

( 

300)  1 

(.000) 

(.000) 

( 

300) 

(•' 

300) 
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TABLE  V 


Same  as  TABLE  III  with  a  total  collision  energy  of  1.75  ev.  Although 
there  are  27  open  channels  at  this  energy,  only  results  for  the  first 
14  are  shown  here.  One  can  see  from  FIG.  V  that  the  probability  of  an 
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FIG.  II 

Intermolecular  Lennard-Jones  potential  for  model  diatom-diatom  system, 
energy  spacings  for  the  collision  partners  and  the  first  several  energy 
levels  for  the  combined  system  are  shown  on  the  same  energy  scale 
for  comparison. 
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Q  -  QUANTUM  MECHANICAL 
S  -  SEMICLASS1CAL 
C  -  CLASSICAL 


0  1  234  56  789  10  1 

INITIAL  QUANTUM  STATE 

FIG.  Ill 

Total  probability  of  an  inelastic  event  for  a  collision  energy  of  1.25  ev. 
Shown  are  probabilities  predicted  by  quantum  mechanical  close  coupling 
calculations  (Q),  seoiclassical  classical  path  calculations  (S)  and 
classical  trajectory  calculations  using  the  histogram  binning  technique  (C) 
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PROBABILITY  OF  TRANSITION  FROM  (1.0) 


01  234  56789  10 

FINAL  QUANTUM  STATE 
FIG.  VI 

State-to-state  transition  probabilities  for  a  1.75  ev  collision 
beginning  in  the  second  entrance  channel. 
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CLASSICAL  ACTION 
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FINAL  CLASSICAL  ACTION  (N-N) 

FIG.  IX 

Distribution  of  final  values  of  classical  action  variables  at  2.25  ev.. 
The  initial  values  of  the  action  variables  for  all  trajectories 
are  oxygen:  2.0,  nitrogen:  1.0  corresponding  to  quantum  state  8. 

The  initial  phases  of  the  oscillators  are  chosen  randomly  between 
0.  and  2  .  The  grid  represents  the  assignment  to  quantum  states. 
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